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Currently there exists no single unified approach for efficiently and accurately solving 
computational fluid dynamics (CFD) problems across the Mach number regime, from truly 
low speed incompressible flows to hypersonic speeds. There are several CFD codes that have 
evolved into sophisticated prediction tools with a wide variety of features including multi- 
block capabilities, generalized chemistry and thermodynamics models among other features. 
However, as these codes evolve, the demand placed on the end user also increases simply 
because of the myriad of features that are incorporated into these codes. In order for a user to 
be able to solve a wide range of problems, several codes may be needed requiring the user to 
be familiar with the intricacies of each code and their rather complicated input files. 
Moreover, the cost of training users and maintaining several codes becomes prohibitive. 


Objective of the Work 


The objective of the current work is to extend the compressible, characteristic-based, 
thermochemical nonequilibrium Navier-Stokes code GASP to veTy low speed flows and 
simultaneously improve convergence at all speeds. Before this work began, the practical 
speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a 
number of new techniques have been developed for more accurate physical and numerical 
modeling. 

A pproach Used 


The primary focus has been on the development of optimal preconditioning techniques for the 
Euler and the Navier-Stokes equations with general finite-rate chemistry models and both 
equilibrium and nonequilibrium thermodynamics models. We began with the work of Van 
Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to 
include three-dimensional reacting flows. The basic steps required to accomplish this task 
were a transformation to steam-aligned coordinates, the formulation of the preconditioning 
matrix, incorporation into both explicit and implicit temporal integration schemes, and 
modification of the numerical flux formulae. In addition, we improved the convergence rate of 
the implicit time integration schemes in GASP through the use of inner iteration strategies 
and the use of the GMRES (General Minimized RESisual) which belongs to the class of 
algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the 
practical utility of GASP through the addition of mesh sequencing, a technique in which 
computations begin on a coarse grid and get interpolated onto successively finer grids. 

Conclusions Relevant to Rocket P ropulsion 


The fluid dynamic problems of interest to the propulsion community involve complex flow 
physics spanning different velocity regimes and possibly involving chemical reactions. This 
class of problems results in widely disparate time scales causing numerical stiffness. Even in 
the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very 
low speed flows which can be quantified by the large condition number of the system and 
evidenced by slow convergence rates. This results in the need . for thorough numerical 
analysis and subsequent implementation of sophisticated numerical techniques for these 
difficult yet practical problems. As a result of this work, we have been able to extend the 
range of applicablity of compressible codes to very low speed inviscid flows (M=.001) and 
reacting flows. Our work now centers on the extension to viscous flows. 
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Constructs syntactically correct input 


Phase I Objectives 
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Transonic Flow over Circular Arc 
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Tangential Injection of Helium into Air 
with Shock/Boundary Layer Interaction 
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Table 1: AeroSoft Phase II Timetable 
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